Strong electron-lattice coupling as the mechanism behind charge 
density wave transformations in transition-metal dichalcogenides 



We consider single band of conduction electrons interacting with displacements of the transitional 
ions. In the classical regime strong enough coupling transforms the harmonic elastic energy for an 
ion to the one of the well with two deep minima, so that the system is described in terms of Ising 
spins. Inter-site interactions order spins at lower temperatures. Extension to the quantum regime 
is discussed. Below the CDW-transition the energy spectrum of electrons remains metallic because 
the structural vector Q and the FS sizes are not related. Large values of the CDW gap seen in 
the tunneling experiments correspond to the energy of the minima in the electron-ion two-well 
complex. The gap is defined through the density of states (DOS) inside the electronic bands below 
the CDW transition. We focus mainly on electronic properties of transition-metal dichalcogenides. 
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I. INTRODUCTION 



Origin of charge density waves (CDW) in the transition-metals dichalcogenides (TMDC) 
is subject of debates since their discovery in the early 70's [1]. At the time, the most 
popular theoretical scenario was the so-called "nesting", the congruency of two or more 
Fermi surfaces (FS) separated by a vector Q in the momentum space. Interactions with the 
momentum transfer Q would lead to the CDW instability (with Q becoming the structural 
vector) and to the opening of energy gaps on FS. 

n 

The alternative explanation [2] ascribes the structural instability to the presence of saddle 
points in the electronic spectrum near the Fermi energy. The logarithmic singularities in the 
electronic density of states (DOS) at the saddle points favor instabilities with the momentum 
transfer, Q connecting the two saddle points. The mechanism also leads to the energy gaps 
in the electron spectrum. 

Subsequent experiments and band structure calculations for real materials gave no sup- 
port to either of the two concepts. The CDW transition does not affect properties of dichalco- 
genides any noticeable. Materials remain metallic below the transition temperature, Tcdw- 
Superconductivity in 2/f-TaSe2, 2/J-TaS2 and 2/7-NbSe2 takes place on the background of 
the CDW phase. The CDW gap values ~ 0.1 eV surprisingly large in comparison with the 
values of transition temperatures, TcDw(Ta) = 122 K and TcDw(Nb) = 33.5 K, were ob- 
served via the dl/dV characteristics in the tunneling experiments jsl. The gap of such order 
of magnitude is noticeable in the infrared data {4] (2if-TaSe2) even close to the transition 
temperature. 

The phenomenological analysis in jsl revealed the short coherence length in the CDW 
phase, ^0 ~ 3 — 10 A, thus implying the important role of fluctuations. 

The density functional calculations (2if-NbSe2) have shown that the contribution from 
the nested pieces of FS into the charge susceptibihty is negligible at the expected value of 



the vector Q, although, indeed, the ground state of the system at T = 0, according to 
is the CDW phase. 

It was well known since 1958 that at the strong enough e-ph coupling the square of 
renormalized phonon frequency turns out negative, thus signaling instability of the lattice 



a. 



In 



through the example of lT-TaS2, it was demonstrated that CDW transitions in 



the transition-metal dichalcogenides can be driven by the Migdal's instability The new 



tinge added to this concept in [8| was that dispersion of the renormahzed (imaginary) soft 
mode can become strongly peaked near the structural vector of instability, Q if the e-ph 
matrix element depends on the transferred momentum. 



There is no sma' 



1 physical parameter in the problem |7[. Correspondingly, there are 



fl. 



no other means in 
numerical calculations. 



8| to discuss electronic properties of 2iJ-TaSe2 and 2if-NbSe2 beside 



As it was already mentioned above, the nesting scenario and that one of Ref. I2I both 



contradict to the experimental data [3|, |J] and to their analysis and the interpretation [5| 
In spite of that, the two concepts are commonly used for interpretation of results even in 
recent ARPES experiments (see, e.g., 

,|iol)- 

We address peculiarities of transition-metal dichalcogenides from the point of view dif- 
fering from jsj. We argue that instead of development of the Migdal instability in the soft 
mode scenario jsj for the propagating phonon mode, the CDW instability may realize itself 
in the two stages. Namely, at first, strong e-ph interactions bring electrons and ions close 



together in kind of a polaronic effect similar to the one first discussed for VsSi in The 
potential in which the ion moves becomes anharmonic. At a strong enough e-ph coupling the 
potential possesses two minima, so that the system can be described in terms of Ising spins. 
The subsequent ordering of the local sites at a lower temperature occurs due to inter-site 
interactions. 



In Ref. lUJ electron-lattice interactions were limited to the interaction of electrons with 
single dispersionless (Einstein) optical mode. In the approach below conduction electrons 
interact with arbitrary ionic displacements. Strength of the intersite interactions and the 
value of contributions into the local elastic potential are separately evaluated. 

"Trapping" of the electronic cloud near an ion by the elastic field is deemed responsible 
for the main energy gains and pre-determines, whereby, by the order of magnitude, the 
energy scale of the CDW gap. According to our estimates, the wells' minima are deep (~ 
few tenths of 1 eV). Thereby, the Kondo-like regime of quantum tunneling that was the 
subject of the main concern in [Hj] seems to be of no relevance to dichalcogenides, at least 
at the temperatures of interest. We focus in the following on properties of the electronic 
sub-system. 



II. CHOICE OF THE MODEL 



Before proceeding further, we briefly summarize the information pertaining to the CDW 
transitions in TMDC 

Symmetry of the order parameters governing the transition was discussed in js], [l2|. The 
lattice superstructure below Tcdw [i| is formed by the triple-Q modulations of the ionic 
positions along the three symmetry axes in the hexagonal 2H-pha.se. The Q-vector may be 
incommensurate, but even then Q remains very close to a*/3 (here a* = A7r/\/3a). 

It was perceived already in the early 70's 13| that the structural changes in 2if-NbSe2 are 



related to displacements of the niobium ions. Correspondingly, in what follows, we choose a 
simplified model and consider only shifts of the transition-metals ions (Nb or Ta) along each 
of these three symmetry lines. The dramatic softening of the longitudinal acoustic mode 



with Q about a*/3 was recently observed in 1^. Note, in passing, that the optical modes 



and the longitudinal acoustic branch of the same Si- symmetry are actually linearly coupled 
near point a*/3 of the Brillouin Zone (BZ). We return to this later to show that softening 
of the optical mode is accompanied by softening of the acoustic phonon, and vice versa. 

Transition-metals dichalcogenides are layered compounds with the weak Van der Waals 
coupling between the layers. We restrict ourselves by the quasi-two-dimensional properties 
(Q2D) of a single layer. 

We consider a single isotropic band of electrons interacting strongly with displacements 
of the transition-metal ions. Implicitly, the band is assumed to be predominantly of the 4f 
(Nb) - or 5f (Ta) - character. 

In the model [8|] the "bare" frequency of the propagating phonon mode, Wo(^) is renormal- 



k ^ Q). Unlike js], our attention is concentrated on the local environment of the single ion. 
In this respect, as we discuss later, our model shows some features common to the Holstein 



ized by the polarization operator 15|] (with the e-ph matrix element g{p, p -|- k) is peaked at 



model 



16|. 



For a heavy ion one may expect that quantum effects are not crucial for properties of the 
lattice, at least, at not too low a temperature. The temperature of the (incommensurate) 
transition in 2iJ-TaSe2, Tcdw = 122 K |l| is rather high and, in the first approximation, 
one can neglect the kinetic energy of the Ta ions. 

With the notations Ui for a displacement of the single ion at the point Ri the "bare" 



elastic matrix is: U{Ri — Rj)uiUj. At i = j it is merely the potential energy of the oscillator: 

U{Ri - R,)uiUj ]^MujIuI = ]^ku\ . (1) 

Interactions with electrons change the elastic matrix. We calculate the contribution to the 
elastic energy that is due to interactions of electrons with arbitrary static displacements. 
With the Hamiltonian for the e-ph interaction in the form 

He-ph = J2gU^^^{RMR^) (2) 

i 

the energy of the coupled electron-lattice system E{ui...un] g) can be calculated from the 
equation: 

dE{ui...UN) \^ /ox /o\ 

— = }^u^n{Ri) , (3) 

where n{Ri) = {'ip^ {Ri)-ip{Ri)) is the number of electrons per unit cell at the point Ri. In 
terms of the electronic Green function (at T = 0, ll5|), Ga/3{R, t; R', t') = —i{T{ipa{Ri, t)t(j^{R', t'))), 



one has: n{Ri) = —2iG{R,t] R,t + 6). (The Green function is diagonal in the spin indices: 

Gai^(^R,t] R' ,t') = — idoif^Gi^Ryt] R' ,t'^ ). 

III. COUPLED ELECTRON-LATTICE SYSTEM IN THE CLASSICAL REGIME 

At static displacements, one needs to know only the frequency component, G{R, R';co): 

— G{R, R'; u) exp{-iuj{t - t')) . (4) 
Re-arrange the power expansion of G{R, R'; u) in Ui as: 

G{R, R'; oo) = Go{R - R'; 00) + ^oiR - Rf, uj)guiGo{R - Rf, 00) (5) 

i 

+ ^Go{R- Ru uj)guiGo{Ri - Rk, uj)gukGo{Rk - R';uj) + ... . 

i^k 

The e-ph contributions in all powers in g are now summed first for the single site: Ui ^ Ui. 
Each of Ui,Uk, ... in Eq. stands for the "dressed" local deformation at the corresponding 
site, Ri, Rk, and is determined by the relation: 

barui =Ui + UiGo{Ri = Rf, u)gui . (6) 



The free Green function, Go{Ri = Ri]u) equals (l5|: 



Go{Ri - Ri] uj) - I ^ — — — -— -: . (7) 



With the help of identity 



P ( ) — msig'n.{uj)5{uj — ^) 



— ^ + 2sign((X')5 \uj — ^ 
and assuming the electron-hole symmetry for |,^| = \E{p) — Ep\, one obtains: 

Go{Ri = Ri] u) = -iTcu{Ep)sign{u) . (8) 

After trivial calculations, one finds the electronic contribution into the local elastic matrix 
at the site Ri, E{ui) = —2i dgUiG{Ri = Ri\ t = t' + S; g): 

E{ui) = -^\n{l + [niy{EF)gu,f} , (9) 

where W is the bandwidth. E{ui) is determined as the energy per one ion, ^{Ep) is the 
number of states at the Fermi level per unit cell: v{Ef) = ) 5*0, where Sq is the area of 
the 2D unit cell. WulEp) = s is an insignificant model parameter, s ~ 1. We take s = 2. 
The interaction between two sites, E{ui,Uk), is given by the second term in Eq. (j5]): 

/did 
— ,Gl{Ri,k;u) (10) 
In 

where Ri^k = \Ri - Rk\- 

The analytic form of Go{Ri^k','^) being cumbersome in 2D, for the estimate we use its 
asymptotic at pFRi,k > 1: 



Go{R;u)) =^ ^'^(-^^H/^T^ ^ isign(a;) 



(11) 



We obtain: 



E{ui,Uk) = UiUkg y [Ef)pfVf—, — 77-5— • (12) 

[PfRY 

The total elastic energy of a single ion is the sum of Eqs. ([1]) and (|9]). Introducing: 

/=^^^V^ (13) 

\2tiu{Ef) ) ^ ' 

- {MuIv{Ef)) = , (14) 



and the dimensionless notations for the ions' shifts, Ui = (ui/uo), one writes down the local 
elastic energy in the following simple form: 

1 



[^?-ln(l + A\^)] 



(15) 



7Tl'{Ep 

Here is the square of the dimensionless e-ph coupling constant and Uq determines the 
spatial scale of the local elastic potential. The energy scale T* = I/ttu^Ef) = W/2tt is 
expressed through DOS in the electronic band. For 2if-NbSe2 the band calculations 
gave v{Ep) = 2.8 states/eV per two bands. Here T* is of the order of tenths of 1 eV. 
At > 1 the potential Utot{ui) has two deep minima at u+- = ±Vl — A"^: 



'1/hv{Ef)) [l-A-2-lnA2] . 



(16) 



At A^ < 1 and temperatures below T* ^ 1 and Utot{ui) can be written as: 



f/tot(«.) = (lAK^F)){^.'(l-A2) + Y^f} . (17) 

The quartic term in Eq. (fT7|) is small, but the anharmonic contribution into the elastic 
ener gy i s the necessary ingredient in the molecular field approach to a CDW transition (see, 
e.g., lisl). 

Re- writing PfVf as: PfVf = {1/txv{Ef)){pISq/2), where S'o = {V^/^)a'^ is the area of 
the triangular unit cell, Eq. (fT2l) can be written down in the notations of Eqs. fll3lll4p : 



E(ui,Uk) 



(18) 



7rz/(£'i7' 

The numerical factor in this expression shows that the inter-site interactions are weak com 
pared to the on-site Utot{ui). 

At A^ > 1 the model reduces to the model of interacting Ising spins: 



u+-{i) =^ a/1 — A 2 cTj, cTj = ±1 . 

When A^ < 1, the expressions for t/tot(wi) and Eq. (ITSl) for E{ui,Uk) together complete 
formulation of the problem: in the classical regime the system is fully described by the 
partition function Z{T,g): 



Z{T,g) = / (Jldui) exp 



(19) 



where Utot{ui,Uk) = Utot{ui) + E{ui,Uk). 

The phase transitions are habitually treated in the molecular field approximation (see, 
e.g., jl8|). The method being not exact for short-ranged interaction; we do not discuss 
details of the CDW transitions itself. 

For local properties one has for the partition function, Z{T,g): 

Z,{T,g) = I dn,exp|-— [u^ - \n{l + A'u^)]'^ . (20) 

Scattering of electrons on the lattice displacements above Tcdw is characterized by the 
average (gui) that enters the denominator of the Green function: 

G-\R; 00) = 00- E{p) + Ef + {gUi) . (21) 

From Eq. ([6]) it follows: gUi = gui[l — isign{uj)TTi'{Ep) gui]^^ . In the normal phase terms 
that are odd in Ui can be omitted. In dimensionless variables: 

1 A^u^ 

gUi = zsign(w) — — — . (22) 

t:v[Ef) 1 + J\. Ui 

At < 1 the imaginary part is: G"^(i?;cj) is {gui) = isign{uj)A'^{n /2)u{Ef)T. In the 
classical regime above Tcdw the resistivity of the system would be linear in T. In the 
opposite limit of A^ > 1, -u^ = _ = 1 — A"^, and the imaginary part is a constant of order 
of 1/ttv{Ef)- The entropy for the system of non- interacting Ising spins is finite. 

In the ordered state with all ions occupying same minima, non-zero gUi 7^ stands 
together with the chemical potential: 

1 



gUi = ±_-— Va^ . (23) 

So far, for simplicity of the arguments, it was tacitly implied that the order parameter 
gUi 7^ in Eq. f l22|) stands for the CDW transition with the structural vector Q = 0. As in 
TMDC the Q-vector is non-zero, the order parameter gUi{Q) in Eq. ( 12T|) couples electronic 
states with the energies E{p) and E{p + Q). 

In such a way, non-zero gUi{Q) changes the energy spectrum and, hence, DOS in the 
vicinity of such paired points, but does it mainly for energies away from the Fermi level 
because in dichalcogenides the Q-vector and sizes of the Fermi-surfaces are not related. Of 
course some Fermi- surfaces points may be affected at onset of the CDW order, as it is really 



observed in js]. 
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IV. DISCUSSION OF THE RESULTS 



A. CDW instability: above and below Tcdw 

It is interesting to discuss the relevance of the above model to the experimental results, 
in particular, for 2/7-TaSe2. 

First, recall that the non-linear potential (|9]) was derived for a single ion. For Eq. (Q 
to be meaningful, uq must be small compared to the lattice parameter, a ^ 3.45 A. With 
the phonon frequencies, Uq for NbSe2 typically in the range 10 — 20 meV, is indeed small: 
Mo ^ 0.30 — 0.15 A. Similar estimates for 2if-TaSe2 give Uq ^ 0.2 — 0.1 A. This justifies the 
assum ptio n. Experimentally, the superlattices shifts of cations are between 0.1 A [l| and 
0.5 A [11 1 and, hence, have same order of magnitude as mq, as it should be in the model of 
Ising spins. 

In y the large CDW gaps ~ 90 meV and ~ 60 meV for 2i7-TaSe2 and 2i7-NbSe2, 
respectively, were derived from the tunneling non-linear current characteristics, dl/dV. Such 
values are in the right correspondence with the energy scale, T* = l/Tiv{Ep) equal to 
0.1 — 0.3 eV, depending on the number of states per unit cell in the specific material (for 
NbSe2 ~ 2.8 eV~^ per unit cell per two Nb bands |17]). So far, however, the CDW gap itself 
was rather vaguely defined. There, it was shown that the CDW with non-zero Q couples 
the electronic states with E{p) and E{p + Q). For the two symmetric points, E{p — Q/2) 
and E{p + Q/2), say, inside the energy band for the Fermi surface centered at the F-point, 
the energy spectrum in their vicinity is: 

Ei,2(p±g/2) = E{p±Q/2) ± \gu,{Q)\ . (24) 

Eq. ( 12^ defines the CDW gap as it appears in the I — V characteristics which measure the 
energy dependence of DOS. The CDW gaps depend on the value of the parameter, > 1 
as in Eqs. (fT5|) and (|22|). 

In 2if-NbSe2 the vector Q ^ a*/3 is shorter than the radius of the Fermi surface centered 
at the F-point. Changes in DOS at the CDW transition for the three symmetric pairs 
of points inside the Fermi surface were directly detected in the ARPES experiment jl^. 



(Actually, in 19| the map of all the spots paired with the Q-vector was obtained for the 
whole BZ; the pattern has the hexagonal symmetry). 

Recall now that theoretically the frequency of the collective excitations behaves differently 



at T = TcDW for the local potentials with one or two minima [20| (see also in [18|]). In the 
former case the frequency vanishes at the temperature of the transition. In the latter, the 
frequency remains finite as an ion is now "trapped" by one of the two minimum. While 
neutron experiments [l](b) for 2iJ-TaSe2 gave finite oo'^iQ) ~ 20 meV^ at Tcdw = 122 K, 
softening of the acoustic phonons was observed for 2iJ-NbSe2 at Tcdw = 33.5 K jl^ . 

The difference in the phonon modes' behavior in the two materials needs a clarification. 
Recall that the longitudinal acoustic branch is coupled linearly with the optical mode, u 
of the same Si- symmetry at Q = a*/3. Free energy then has a contribution of the form: 
F{u,s) = Ug{s'^/2) + a;^(n^/2) + tus, where s stands for the acoustic branch. Minimizing 
F{u,s) and excluding u gives F{u,s) = [w^ — (t^/u;^)](s^/2). Let the optical mode, u be 
the mode that drives the transition. At ^^(T) (one minimum), the acoustic mode is 
the first one that manifests the onset of the transition. If the potential has a few minima, 
the effective frequency of the acoustic mode oj^^q = [w^ — (t^/w^)] at Tcdw may or may not 
be zero depending on the temperature behavior of ujI{T). Data Ll](b) for 2if-TaSe2 agree 
better with the second possibility. 



The two 2/J— materials have different masses of Ta- and Nb- ions. It is known [20] that 
at lower temperatures, when quantum effects prevail, ^^(T) — )■ even for a potential with a 
few minima. This could be another possible interpretation for the acoustic phonon frequency 
vanishing at Tcdw = 33.5 K in the 2if-NbSe2 1^ and the finite ui'^{Q) at Tcdw = 122 K in 
2i7-TaSe2 Q(b). 

The thermodynamics of the second order phase transition would look much alike for 
the displacive transition or for the ordering of the Ising spins. The finite uj^{Q) at 
Tcdw = 122 K in 2if-TaSe2 [l|(b) is the first argument supporting the Ising spins model. 

We argue that the temperature dependence of resistivity in all TMDC's seems to indicate 
in the same direction. In fact, for a broad temperature interval the resistivity, p(T) behaves 



as p(T) = po + clT, with the large intercept po ~ 100 — 150 pi7 • cm (2if-TaSe2, |2l|(a. 
b)). This behavior is consistent with the inverse mean free time in Eq. fl^T]) for the Green 
function of the form r^^ = Tq ^ + dT and = {1/tti'{Ef)){1 — A^^). 

In the ARPES experiment [lO( the self-energy, Ti{u) = S'(a;) + iT,"{u) was measured 
directly on the Fermi surface of 2if-TaSe2 centered at the F-point of BZ. The imaginary 
part, S"((X') ~ 60 — 70 meV is frequency independent above Tcdw; this value is consistent 
with Tq^ T* = l/niy^Ep). An estimate with Tq^ ~ T* gives the right value of po from 
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2l|(a,b). 

As to the real part, is negligible at T = 111 K (i.e. above Tcdw = 88 K). 

Below Tcdw ^'{^) starts to show the behavior typical for the self-energy of normal electrons 
interacting with phonons [7] and is peaked at the phonons frequencies ~ 50 meV. Such typical 
metallic signature obviously agrees with the suggestion that the CDW order parameter 
gUi{Q) in Eq. fl23|) affects the electronic states of the energy band at the F-point only below 
the Fermi level 221. |23|. 

Finally, the NMR methods may help with revealing properties of the transition atoms 
occupying the symmetry positions in the commensurate phase of 2if-TaSe2, as it is pointed 
out below. 



B. Quantum regime 

At low T quantum effects become important in few aspects. Consider first the Schrodinger 
equation for an ion moving in the rigid potential f/tot(^j) = {^/T<'^{Ep))U{ui), with U{ui) in 
the dimensionless notations. With the kinetic energy in the same notations, one has: 

^ -^<il{u) + [U{u)-E]^{u) = (25) 



2M du"^ 

The dimensionless "mass", M in fl25|) is defined by 1/M = ['Kh'{Ep)uJo]'^ /2. (The adiabatic 
parameter 1/y/M = {tiv{Ep)ujq) / \/2 is about one tenth in 2if-NbSe2). 

For the single-minimum well (A^ < 1) quantum effects are important when Tcdw < ^^o- 
The lattice oscillations are now quantized and the electronic Green function is "dressed" by 
phonons. In the adiabatic approximations, solution for the problem of interacting electrons 
and phonons in normal metals was given many years ago in [3]. 

At first glance. Eg. ( 125|) adds a possibility of quantum tunneling between minima of the 



two-well potential [20|. Actually, Eq. (125|1 could account only for tunneling in a "rigid" 
potential, a potential built in the lattice only by the inter-atomic forces. The two-well 
potential (Eqs. ( I91IT51) ) is formed by interactions between electrons and the local lattice 



distortions. It was emphasized in [HI] in connection with the martensitic transitions in 
NbaSn and VsSi, that at the tunneling event, at which the ion goes over, say, from 
to M_, the electronic configuration reverses as well. Such a feature cannot be described in 



terms of Eq. (I25l) . According to ll|], the under-barrier tunneling results in a Kondo-like 
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quantum regime at which the height of the barrier would diminish with temperature. The 
parameter Vm stands in the exponent of the expression for the tunnehng matrix element. 



With no closed solution found in one may only argue that with the value of v M ~ 10 
the Kondo-like regime [ll| is probably of no relevance to the CDW physics in TMDC 

Below the transition in the CDW ground state ions will occupy their proper minima, in 
accordance with the superlattice pattern [l!](b) and tunneling between minima must stop. 
Nevertheless, experimentally, the pattern reveals an interesting peculiarity. Indeed, in the 
CDW phase one in three atoms along the symmetrjy lines finds itself in the position with the 
trigonal symmetry (2i7-TaSe2 [li(b); see Fig. 4 in [12]). While below transition the intersite 
interactions do indeed arrest quantum tunneling between minima for the other two of the 
three atoms, the degeneracy is not lifted for the atom in this symmetric position. 

NMR experiments seem be able to verify the very concept of the two minima-potential 
by studying the . positions with the tngonal sy..et.y 

According to ARPES data [9], the CDW transition in 2iJ-NbSe2 slightly affects only 
the two-barrel Fermi surface at the K-point in BZ, while FS at the F-point remains intact. 
The three Q-vectors couple together the three points on the inner FS at K, which are seen 
experimentally with small but observable gaps at these points. Surprisingly, the observation 
of such small local gaps (~ 2.4 meV) was interpreted in {q] as the conclusive proof in favor 
of the nesting mechanism of the CDW formation. 



The physics of strong local e-ph interaction bring us back to the Holstein model [16| of 
electrons interacting with dispersionless phonons. No exact solution is known for the Hol- 
stein model either, but its low teraperature physics was investigated numerically in DMFT 



(Dynamical Mean Field Theory, 2m) approximation. (For a brief summary of results for 
the Holstein model at T = see j25[)- 

The DMFT approach is strictly local. Intersite interactions (see Eq. (fT8l) ) and the CDW 
transition itself cannot be treated by DMFT. Here we indicate parallels between our physics 
above and the results at T = for the Holstein model. Among them are: (1) The 
double-well potential that develops when the e-ph coupling constant becomes larger some 
critical value; (2) Large imaginary part in the Green function: Spectral Function extends 
over the energy interval that significantly exceeds the phonon frequency (compare with our 
Eq. ( 122])) ; (3) The ground state remains metallic (at least for not-too-strong e-ph coupling). 

In 251] the mass of electronic excitations increases as the residue at the pole of the 
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electronic Green function decreases. 

Judging by these results, the physics studied above allows its extension into the low 
temperature regime. 

V. SUMMARY 

In summary, we calculated the local elastic energy for single ion by re-summing exactly 
the electron-lattice interactions in the real space. Strong e-ph interactions qualitatively 
change the local environment by binding ions and electrons together, thus breaking the 
adiabatic approximation. The concept is the realization of the Migdal instability different 
from that one considered in [8,]. The CDW transition takes the form of a phase transition in 
the system of interacting Ising spins. The value of the structural vector Q and parameters 
of the Fermi surfaces being not related with each other, the energy spectrum of electrons 
remains metallic below the temperature of the CDW transition. The CDW gap seen in the 
tunneling experiments below the transition is defined in terms of reduced DOS inside the 
electronic bands. Its large value is consistent with the energies of the deep minima of the local 
potential. Other experimental evidences in favor of the Ising model were enumerated. The 
electronic properties of the transition- metal dichalcogenides, as they are seen by ARPES, 
agree well with the suggested concept. 

Conclusions from the analytical results derived in the classical regime can be extended to 
lower temperatures as it follows from the comparison with numerical results for the Holstein 
model. 
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